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ABSTRACT 

The set of moments and the distribution function of the one-electron current in a one- 
dimensional disordered ring with arbitrary magnetic flux are calculated. 



★ Also Sezione INFN di Milano, 20133 Milano, Italy. E.mail: cattaneo@vaxmi.mi.infn.it. 
f Also Sezione INFN di Pavia, 27100 Pavia, Italy. E.mail: gamba@polito.it. 

| Permanent address: Budker Institute of Nuclear Physics, 630090 Novosibirsk, Russia. E.mail 
kolokolov@vaxmi. mi. infn. it 



1. Introduction 



The theoretical work of Ref. 1, together with recent experimental results [2, 3] have stirred 
a strong interest toward the problem of the persistent current in a mesoscopic metal ring 
immersed in a magnetic field. This current has been computed using various approaches 
[4-11]. In some works [4-7] it was considered the experimental situation of a ring of finite 
thickness, such that the number of transverse channels be much greater than one. This 
allowed the authors to use the methods of weak-localization theory. In this case one has also 
to take into account the electron-electron interaction [8]. 

On the other hand, the computation of the one-electron current in an idealized one- 
dimensional disordered ring is also of interest, at least from the theoretical point of view. 
One-dimensional localization effects lead in this case to a non-trivial current dependence on 
the magnetic flux (see Ref. 9 and below). In Ref. 10 such a calculation was performed, 
but only in the weak magnetic field case. In Ref. 11 the one-electron current averaged over 
an ensemble of rings was derived non-perturbatively using Grassmann matrix integration. 
However, the unexplicit form of the resulting expression and the need of tedious computations 
do not allow to check the validity of some approximations. 

In Ref. 9 it was developed a new path integral approach to the study of one-dimensional 
localization. Along with the multipoint density correlators, this new method allowed to 
compute the averaged absolute values (|/|) of the one-electron current / in a disordered 
metal ring with arbitrary magnetic flux. 

In the present paper we show that the method introduced in Ref. 9 allows us to recon- 
struct completely the distribution function V$(I) over an ensemble of one-dimensional rings 
with given magnetic flux $. Simple explicit expressions for the moments (I 271 ') of the current 
/ are also obtained 
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2. Path Integral Representation 

Let us recall the main steps of the path-integral approach introduced in Ref. 9. The 
Schrodinger equation 

(H - k 2 )tp = (-^ + U(x) - k 2 ^ = (2.1) 

maps the two-dimensional space of the initial conditions (tp'(xo) +iki/)(xo), i^'{xq) — ikip(xo)) 
to the two-dimensional space of the solutions at the point x through the matrix 

T(x, x ) = e ik ( x - Xo ^ z T(x, x )e ik ( x - Xo ^ z , (2.2) 

where T(x,xq) obeys the following first-order equation: 

= {itp{x)s z + (+(x)s- + C(x)s + )T, (2.3) 

and 

<p(x) = ~U(x), fyx) = ±^U(x)e ±2 * kx . (2.4) 

Here s z = a z /2 and s ± = (a x ± i<r»)/2 are the usual spin operators. It has been shown [12] 
(see also Ref. 9 for greater detail) that in the limit 

U > 1, (2.5) 

where i is the mean free path, the fields ip(x) and C ( x ) can be considered as statistically 
independent. If the initial potential U(x) is a Gaussian random function of x with correlator 

(U(x)U(x')) = D8(x-x'), (2.6) 

then the averaging weight over the fields (p(x) and C ± ( a; ) nas the form: 

V<p(x) Vfyx) exp (-t ■ J dx l^ 2 (x) + C+(x)C(x)^j , (2.7) 

where £ = Ak 2 /D is the localization length. It has been shown [9,13-15] that the following 



change of variables in (2.7) 

i<p = ip + 2ip + ip~ , 

C = "0" - i(*l>~ - ^ + (^") 2 , 

c+ = ^+ 

brings the operator T(x, xq) in the form of a product of usual matrix exponentials: 



(2.8) 





X 




T(x, xq) = exp [s + tp (#)] exp 


is z J dtp(t) 


exp 









Xo 



(2.9) 



This statement can be checked by deriving the evolution equation for the operator (2.9) and 
comparing it with (2.3). The field ip~(x) is assumed to obey the initial condition tp~(xo) = 0, 
thus providing the equality T{xq,xq) = 1. Under a proper regularization, which is required 
by physical considerations (see Ref. 9), the Jacobian of the transformation (2.8) is seen to 
be equal to 

/ L \ 

(2.10) 



J oc exp I -- / dtp(t) | . 



The surface of integration in the space of the complex fields (p, ^ is defined by the equations 

lmp = 0, C = (C + )*, (2.H) 

where * denotes complex conjugation. This surface can be deformed to the standard one 

lmp = 0, ^ = (^ + )*, (2.12) 

if all the quantities to be averaged are written in a form which allows analytical continuation 
from the surface (2.11) to the whole complex space of field configurations (for more details 
see Ref. 9 and 13). This requirement turns out to be fairly constructive. 
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3. Calculation of the Current Moments 
and of the Current Distribution Function 

In an appropriate gauge the wave-function of an electron moving in a metal ring of size 2L 
immersed in a magnetic flux <3>, measured in units of flux quanta, obeys eqn. (2.1). Topology 
and flux dependence are then encoded in the boundary conditions 

(ip'(L) ± ikip(L)) = e 2 " $ (V/(-L) ± ik^{-L)). (3.1) 

The mean value of a function /(/) of the current / can be defined as follows: 

(/(')> = ^E^-^)/(^))> where jn = ~hit- (3 - 2) 

Here E = k 2 is the electron energy and E n are the eigenvalues of the Hamiltonian (2.1) with 
boundary conditions (3.1), which can be written in terms of the matrix T = T(L, — L): 

det(T-e 2 ™ $ ) = 0. (3.3) 

The matrix T = T(L, —L) satisfies the "unitarity" conditions: 

o z t\<j z = T" 1 , detT = l, (3.4) 

and therefore admits the following parametrization: 

/ e* Qs cosh T e lf3a sinh T \ 
T=[ R _ . (3.5) 

V e iPs sinh T e * Qs cosh T J 

Here a S) (3 s and V are by construction (see Ref. 3) slowly varying real functions of L. Sub- 
stituting the parametrization (3.5) into (2.2), we obtain from (3.3) the equation determining 
the set of the eigenvalues E n [10]: 

t(E) = coshT cos(« s + kL) = cos27r$. (3.6) 

Let us start the computation of (/ 2n ): 

. o„. /27rA; v ^ t . / ^ ^ . 9rn \ /2irk r/ , _ _ sin 2m 27r$ \ . „, 

< J > = ( T E^ " E n)jl m ) = ( T W - ^^ p^yppi ) • (3-7) 

Here (3.3) and (3.6) have been taken into account. The (^-function can be eliminated in (3.7) 
using the following consideration: for kL 3> 1 the result of the average (3.7) does not change 



when L varies on a scale much less than t. Then (/ 2m ) must coincide with its average over 
an interval AL of lengths L: 

L+AL 

(l 2m ) L = -£l j dL < /2m >L > where ^ < AL < /. (3.8) 



L 



We can interchange the order of the two averages; then, using the (approximate) constancy 
of the variables T, a s and f3 s on the interval AL, we obtain: 

// 2 ™\ = (J sin 27r$) 2m / = — \ , (3.9) 

V 1 l ° ^ \(sinh 2 r + sin 2 27r$)^/' 1 ' 

where we have set Io = 2k/ L. Eqn. (3.9) can be rewritten in the form: 

00 

(I 2m ) = sin 27r$) 2m J dnn 2 ™- 1 ^-Asmh^r+sin^^ (3 1Q) 



It is important to notice that sinh T can be expressed in terms of the elements of the matrix 
T without using any complex conjugation: 

sinh 2 r = (l 0)7VtQ^, (3.11) 

where t denotes the usual matrix transposition. Thus the above mentioned analytic contin- 
uation from the surface (2.11) is possible. Substituting (2.9) into (3.10) we obtain: 

L 

sinh 2 r = V~(L) J dt^^e-^t^'p^'). (3.12) 

—L 

The average (3.10) is performed using the weight 

VpV^ e- s ' { P^\ (3.13) 

where the action S\p,il) ± ) is obtained from (2.7) after the substitution (2.8), taking into 



#1 Such a procedure, which was proposed in Ref. 9, seems to be equivalent to the rings-ensemble averaging 
of Ref. 11. See also Ref. 16. 



account the Jacobian (2.10), and reads 



S'ip,^) =1 J dx 

-L 



1-3 3 

-p 2 + ^ + ip~ - -ipip+vp- - -(ip + ip~) 2 



-L 



+ - I dxp. (3.14) 



This action has the form of a (0 + 1) -dimensional Schwinger model and the D^-integration 
in (3.10) can be performed using the so-called "bosonization" method [17], representing 



cxp 



dx 



—L 



Vr] exp 



lj 

~e J dx (r] 2 + 2#+^-) 



—L 



(3.15) 



and: 



e -M 2 sinh^r = 1 I dzdz * e -\z? 
7T 



—ipztl) (L) — ipz 



Lj 

J dxi[) + \ 



x)e 



-iJ x L dt P (t) 



(3.16) 



We eliminate the pip + ip and rjip + ip interaction terms through the following gauge trans- 
formation: 



ip{x) — > ip (x) exp 



X 

±1 J dt (277 -ip) 



(3.17) 



which has the Jacobian 



Jr oc exp 



~J dt(2r]-ip) 



(3.18) 



The -integration becomes Gaussian and can be easily performed. Introducing the vari- 
able £(x) and denoting the x-derivative with a dot we have 



f = -3r/ + -p, £(L) = 0, Vp V V oc Vp Vt, 



(3.19) 
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and performing the Gaussian Dp-integration, we come to the following expression for (/ 2m ): 



oo 

(I 2m ) = 2 (J S in 2vr$) 2m [ dfifi 2 ™- 1 [ dzdz 
7r(m-l)! J J 



Lj 

J Vi exp -- J dx(i 2 + ^ 2 \z\ 2 e-t) 



* -fi 2 sin 2 27r<E>— 1^| 2 _ 



■Me ™ J £>£ exp 

t{L)=Q 

£2m 



-L 



e 2 



(J sin 2tt$) e 



in terms of usual quantum-mechanical matrix elements with Hamiltonian 

1 d 2 



H = - 

id? 4 



(3.20) 



(3.21) 



where the ket- and bra- wave funtions are: 

T(0 = e-§, lf%r) = e^(— V2)exp p 2 ^*" ) 



(3.22) 



The factor A/" in (3.20) is the normalization of the standard Feynman-Kac path-integral; 
together with e~27 it provides the equality (1) = 1. Using the complete set of eigenfunctions 
of H 

/,(0 = fVz/sinh27n/i^ [le~i) , tf/,(0 = 4 ^(0, (MM = ^ " 

(3.23) 

where is the standard notation for the modified Bessel function, we obtain, after some 
arithmetic: 

_l)m-l f q 



(I 2m ) = (J sin 27r$) 



2m 1 



(m- 1)! V<9(sin 2 27r$) 



m— 1 



J'Vo sin 2tt$)~ 



(3.24) 



and 



< /2 > = 7=7777- 7^3/9 ( 7 sin 2tt$) 2 / 



v^F(2L/£) 3 / 2 
log 



--^x 2 

xe il 



o 



a/ sinh 2 x + sin 2 27r$ 



sinh 2 x ( sinh 2 x sinh 4 x 

1 + 2 , „ _ _ + 2 ( ; . _ + 



1/2' 



sin 2 27r$ \ sin 2 27T$ sin 4 27T$ 



(3.25) 
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In order to reconstruct the distribution function V$(I) we use the following identity: 

(3-26) 



V*{J) = i hm Im JL lim Im f jq 

n=0 v 7 



It is convenient to compute the sum in (3.26) considering it as the result of an analytical 
continuation in J from the imaginary positive semiaxis. Substituting (3.24) into (3.26) we 
see that for such J the summation gives a well-defined translation operator in the variable 
sin 2 27r$ acting on (^I 2 (Iq sin 2ir$)~ 2 ^. Performing the analytical continuation to the real 
axis we obtain 



V*(I) = 0, for |/| > J 



and 
where 



2e a? (J sin 2n$) : 

0F(2L/£) 3 /2 f/p 



+00 



dx 



xe 2i x 



A(7) 



; 

J sinh 2 x — sinh 2 X(I) 



(3.27) 



for \I\ < J , (3.28) 



2k 



Iq = — and X(I) = sinh 1 

Lj 



Sm27T<I>| I y| - 1 



1/2- 



(3.29) 



It can be checked that this distribution function reproduces all the moments (/ 2m ) as well 
as the result for (|/|) obtained in Ref. 9. 
In the limit /^Owe get 

| Iq sin 27r$| 



A ~ log 2 



and 



Aexp(-^(A-^)r(4, + i) 



(3.30) 



(3.31) 



Um/i)^ r(^ + i) 

in the proximity of the maximum A ~ j. In the limit Lji — > +oo we get for the quantity A 



the normal distribution 



r$(i) di 



e 2L 



) 2 



dX. 



(3.32) 



v/2717^ 

Thus in the thermodynamic limit the fluctuations of A are suppressed and A becomes a 
non-random quantity. This fact is in a deep connection with an earlier result [18] about the 



asymptotically normal distribution of the logarithm of the static resistivity (see also Ref. 12). 
In both cases we are dealing with the response of the system to an external field. In our 
case sin 2 27r$ can formally assume an arbitrary value and in some sense we are considering 
a non-linear response. However we see from (3.32) that in the limit Ljl — > +00 the response 
becomes effectively linear. 



It is worth noting that in the limit 27r$ — > 0, when £/L < 00 is fixed, all the moments of 
the current J tend to zero. This seems natural since in a given potential without symmetries 
(£/L < 00) and with zero magnetic field all the stationary states of the electron in the ring can 
be described by real wave functions. The corresponding quantum-mechanical expectation 
values of the current operator are equal to zero. On the other hand, if we take simultaneously 
the limit £/L — > 00 (free motion case) we can obtain a non-zero result. 

Let us also note that the formal substitution Iq — 1, sin 2tt$ = 1 and J 2 -> T in (3.28) 
gives us the distribution function for the transmission coefficient T. In the limit Ljt — > 00 
this reproduces the known results [19,20] for the moments (T n ), but our formula is valid 
for finite values of L/i as well (the only limitation is that the sample length 2L and the 
localization length £ be great in comparison to the wavelength l/k). 



We would like to thank M. Martellini for his interest in this work and for his continuos 
support. We are also grateful to the referee for his observations which helped us to improve 
the text of this paper. 



When J -> J (3.28) gives 
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